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Abstract—The Geostationary Interferometric Microwave 
Sounder (GIMS) is a new concept of atmospheric microwave 
sounder for China’s future geostationary earth orbit 
meteorological satellite (FY-4). A novel self-calibration method 
for interferometric radiometers with rotating thinned array, 
especially for GIMS, has been proposed in this paper. Compared 
with traditional relative phase calibration approach, neither 
dedicated hardware nor dedicated calibration working model is 
needed to achieved the relative phase calibration in this novel 
inter-element phase calibration. The self-calibration approach is 
inherently merged with the nominal observation working model 
of GIMS, thanks to the continuous array rotating of GIMS 
instrument. A running average scheme has been introduced into 
the self-calibration approach to enhance the signal to noise ratio 
(SNR) of the calibration data, which is normally very low with the 
natural earth scene. The method is demonstrated by both 
simulation and field imaging experiment. 


Index Terms—inter-element phase calibration, Geostationary 
Interferometric Microwave Sounder (GIMS), phase error, 
synthetic aperture radiometer, 


I. INTRODUCTION 


Realizing microwave atmosphere sounding from 
geostationary earth orbit (GEO) is an emerging researching 
field since 1990’s'""""!. GEO can guarantee the fixed continuous 
observation to a particular region, while the relative low 
frequency band (compared with infrared waves) can guarantee 
all-weather observation with a 3D image capability for the 
atmospheric temperature and moisture. Those features make 
the GEO microwave sounding very appealing for short term 
meteorological forecasting. 

The largest obstacle of GEO microwave observation is the 
restriction of spatial resolution. Interferometric aperture 
synthesis is an effective solution to this problem, by using a 
thinned antenna array to replace the single large aperture 
antenna. Since 2000’s, several interferometry based instrument 
concepts have been proposed, including the GeoSTAR"! 
concept by NASA/JPL, the GAS"! concept by ESA and the 
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GIMS"! concept by NSSC. 

GIMS is a millimeter wave imaging sounder concept 
proposed for China’s next generation geostationary 
meteorological satellite (FY-4M). GIMS is supposed to utilize 
a rotating circular thinned array, instead of a stationary Y-shape 
array, to reduce the required number of the antenna elements. A 
proof-of-concept 53GHz ground-based GIMS demonstrator 
with 28 elements has been successfully developed and tested 
during 2009~2012'°!. A new GIMS demonstrator development 
plan has been approved for further risk mitigation. A 70 
element demonstrator is foreseen to be finished before the end 
of 2016. Fig. 1 illustrates the artistic view of GIMS instrument 
in orbit. 


Fig. 1 (a) Conceptual view of GIMS 

The calibration of GIMS is a key process to ensure 
quantitative conversion of raw measurement data into scientific 
data with physical meaning. For GIMS, the periodical 
single-element total power calibration can be easily achieved 
by acold sky mirror and a blackbody, which are installed on the 
fixed part of platform with a deployable mechanism and can be 
periodically switched to block the FOV of each element 
antenna”, Besides the single-element total power calibration, 
periodical inter-element phase calibration is also needed for an 
interferometric radiometer like GIMS. The traditional method 
is to use correlated noise injection”, which would seriously 
increase hardware requirements and system complexity", 
Different from this, some phase self-calibration methods have 
been proposed for GeoSTAR®!"""!” in which, specific extra 
external source is necessary to improve measurement signal to 
noise ratio (SNR) in calibration”), 

In this paper, two new schemes have been introduced in 

GIMS inter-element phase calibration method. 

(1). Thanks to the array rotating, each baseline can get a 
conjugated measurement after half array revolution. 
Taken advantage of this feature, running average 
scheme is proposed to reduce the effect of noise in 
phase calibration for GIMS, so extra point source is no 
longer needed. 


(2). Taken advantage of the prior information of the 
brightness temperature distribution as viewed from 
geostationary orbit, the prior knowledge of the very 
short baseline phases is used to eliminate phase 
ambiguity. 

With above improvements, a complete practical 
inter-element phase calibration method has been set-up for 
GIMS. The method is demonstrated by both simulation and 
filed imaging experiment. 


II. GIMS INTER-ELEMENT PHASE MODEL 


As a microwave interferometric radiometer (MIR), the 
thinned antenna array of GIMS with 70 elements is shown in 
Fig. 2(a). Although the snapshot uv samples are very sparse 
and irregular, as shown in Fig. 2(b), the full sampling coverage 
of spatial frequency domain will be achieved after half 
revolution, which is also necessary for inter-element phase 
calibration. 


Aray Layout Snap-shot uv sampling 
— - 


4 3 2 a 0 1 2 3 4 
u(m) 


Fig. 2. (a) Antenna array of GIMS (b) Corresponding snapshot 
uv samples 

The longest baseline of the array in Fig. 2 is about 6177, 

while the maximum spacing between two adjacent antennas is 


controlled under 130 Å , which is an important coefficient for 
inter-element phase calibration. Similar with the 28-element 
array of the proof-of-concept GIMS demonstrator", there are 
two kinds of antennas in the 70-element GIMS array: three 
small aperture antennas to form the shortest baselines, and 67 
large aperture antennas to form other baselines. It should be 
point out that, the second shortest antenna spacing is 6.37, 
which is short enough to be preset with a confidence interval in 


the following sections. 
The output of receivers i and j is given by!!! 
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In which, u and v are the projections over the X-Y axes of 
the vector defined between the two phase centers of the 
antennas normalized to the wavelength. Subscript i and j 


means channel number. (€,77) are the director cosines!!!) 


which respect to the (X,Y) 
(sin@cos@,sin@sin9) . T, is the 
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axes, equal to 


receiver’s physical 


temperature assumed equal in all elements. The 


contribution of T, is negligible in the GIMS prototype due to 


i 
the low antenna coupling and the large antenna separation. 
F(é,1ņ) describes the plural antenna pattern and F, presents 


the characteristics of interferometer receivers, which is known 
X F 5 13 
as fringe washing function". 


Measured visibility phase of V; is given by 


ġ =Y; ta -a +A +a, (2) 
Where y ij is the ideal phase of the visibility. œ, is a random 
phase term introduced by thermal noise. Œ; and @ j are the 


receivers phases. Œ; is the inseparable phase error introduced 


by fringe washing, which can be considered negligible for the 
short and medium length baselines of GIMS. By running 
average scheme, which will be described in the following 


sections, @,, can also be reduced negligible. Then Eq. (2) can 


be further simplified as 


9; =W; t-a; 6) 


HI. INTER-ELEMENT PHASE CALIBRATION 
Lots of redundant information can be achieved during 
observation for GIMS and element phase information will be 
implied in them. Inter-element phase calibration can extract this 
information by redundant observation in a reliable way. 


A. Optimized Redundancy Equations 


Image of the physical world has the property of real value 
when looking from space. Therefore, its ideal Fourier 
transformation is a Hermitian type function that has the 


property 


V (u,v) =V*(-u,—v) (4) 
In which, asterisk means Hermitian operator. Image will be 
assumed changeless during one period of array rotation, then 
V(u,v) and V(—u,—v) are both measured by the same 
baseline with 180° baseline rotation. Submitting the phase 
information of (4) into (3), we can get” 
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f; =W t-a; (5) 

p7 =-yi+a,-a, 
The superscript 6 means baseline rotation angle. M and 
p” are the phase of vs and y/*7 respectively. With a simple 


algebraic transformation 
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However, this procedure would yield large errors in the 
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estimation of phases because the phase Vii from h; and É; 
presents a large error in most visibility samples due to very poor 
measurement SNR. 


Running average scheme is proposed to optimize Eq. (6). 


Selecting an arbitrary time slice during the array rotation, we 


can get many pairs of p and p” . By averaging them 
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where the superscript bracket means the optimization of y, . 
ij 
[9,,9,] is the rotation angle corresponding to the selected time 
slice. Eq. (7) implies that (g; + or) can be regarded as a 
0. ; i 

measurement of Vij? and W; is the weight of this measure. 
multiple measurements will be weighted averaged to improve 
measurement SNR, which can be defined as the radio of 
visibility amplitude to the standard variation of noise. 

From another point of view, running average scheme is 


equivalent to extending integration time, and the standard 
variation of noise is inversely proportional to the square of the 
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integration time. So W; can be optimized under the SNR 
optimization principle as 

2 

w =( ) (8) 

With running average scheme, the measurement SNR will 

not be restricted by measurement amplitude, but the extent of 

selected time slice, which is limited by the instrument stability. 

Taking all of the baselines into account for an array which 


has N channels similarly as Eq. (7), the following set of 
equations can be established 
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This can be written in a simplified form as 
r.=P-A (10) 


When the first channel is set as a reference by assigning a 
zero phase to it. Eq. (10) is proved a determined and complete 
system of equations for solving the unknown phases in a 
least-square sense as long as the equations related to all 
adjacent antenna pairs are preserved”. So the equations 
corresponding to long baselines should be removed because of 
the serious fringe washing or poor SNR. 

In addition, equal length baselines 
redundancy equations during array rotation 
described as 


Va —Vy=(G, -@,)-@,-@) = 65, -O (11) 


could produce 
[10] which can be 


where baseline ab and baseline ij have the same length. 


Similar as Eq. (7), running average scheme and long baselines 
removing are also operated for Eq. (11), the optimized results 
of all equal length baselines can be written in a simplified form 
as 


Î,=P -A (12) 
Eq. (12) is welcome to be added to Eq. (10), so the final 
redundancy equation is achieved: 
Ô| TP 
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B. Phase Unwrapping 


Phase ambiguity is produced by the division operation in Eq. 
(7), so the solution of the final redundancy equation is 


A +kz , in which K is a vector consist of 0 and 1, other 


information must be added to get accurate A. 

For the very short baselines, which normally means the 
shortest and the second shortest baseline, the visibility is 
mainly determined by earth/sky and land/sea_ contours, 
insensitivity to the brightness temperature details of the earth 
disk. So the phase can be predicted in a confidence interval. 
Thus, the phase ambiguity of these baselines can be easy solved 
by a comparison between the solution of final redundancy 
equations and the preset confidence interval. 

For other baselines, some techniques should be considered 
for phase unwrapping. One is closure phase!'*!, which is widely 
used in optical synthetic aperture field to eliminate the phase 
error of atmospheric disturbance. For microwave synthetic 
aperture radiometer, the closure phase relationship can be 
derived from 
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superscript means measurement moment. Baseline ab ,bc 
and ca stick to each other in a closed loop, which is so-called 
phase closure. Eq. (14) implies that, if two of the phase closure 
have been phase unwrapped, the remaining one will be solved 
by closure phase relationship. 

Another technique is equal length baselines, which has 
already involved in final redundancy equations and now is used 
again for phase unwrapping. 


(14) 
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For equal length baselines ab and ij , we can get 


Vin = Vij (15) 

In which, as one of these equal length baselines has been 
phase unwrapped, the remaining baselines will be solved. 

Eq. (14) and Eq. (15) would also yield large errors due to 
very poor measurement SNR, which can be optimized by 
running average scheme similarly as Eq. (7). 

What should be noticed that phase unwrapping can only 
deduce phase ambiguity and have no influence on the solution 
accuracy, which is decided by the final redundancy equations. 


IV. SIMULATION AND EXPERIMENT 


A. Geo-image Simulation 


In order to assist the implementation of GIMS, a simulator 
is developed to calculate the visibility function and 


reconstructed image!'”!, Pseudo-polar Fast Fourier Transform 
Algorithm is used in the simulator!'®! . 

The original input of simulation is a full earth imaging in 
50.3GHz, which is shown in Fig.3(a). The antennas array is 
shown in Fig.2(a). After brightness to visibility simulation, the 
maximum visibility amplitudes for each baseline are shown in 
Fig.3(b). In which, the visibility amplitude drops dramatically 
as the increase of baseline length. Inter-element phase 
calibration has been conducted with the baselines less than 
130A . The corresponding visibility amplitude is about 0.25K. 


as 309 
re 3 
on 
205] 250 
eoo 
ons| 200 d 
o 
as o 048 g 005 a tal i 100 200 300 400 500 600 700 


5 baseline length wavelength 
(a) (b) 
Fig. 3 (a) original input imaging, (b) visibility amplitude 
corresponding to different baselines 
During inversion simulation, thermal noise with standard 
deviation øo =0.06K corresponding to 500ms integration time 
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is add to the visibility. Random phase errors within +40° are 
added to the receiving elements. The reconstructed image 
without phase calibration is shown in Fig.4(a). Annular fuzzy is 
serious because of the phase error. To calibrate this error, 
Snapshot SNR in Eq. (6) and Eq. (11) is just about 4. After 
optimized by running average scheme, in which the averaging 
time slice is selected as 5min. SNR can be improved larger than 
26, which gives a random phase q@, with a standard variation 


lower than o, =1.3°. For the shortest baselines with a length 


3.14 and the second shortest baselines with a length 6.34 , the 
phase will be preset in [-z/2~7m/2] and [7/2 ~3z/2] 
respectively. The reconstructed image after inter-element phase 
calibration is shown in Fig.4(b). 
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Fig. 4(a) raw retrieved imaging, (b) retrieved imaging after 
inter-element phase calibration 
To evaluate the influence of residual element phase error, 
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radiometric error O% is defined as the standard deviation of 


AT(€,,,7,,) > Which is the difference between the ideal and the 
error contained retrieved image inside the circle with radius 
r=0.1, corresponding to the 45° incidence angle from 


GEO"!, Element phase error O p is defined as the standard 


deviation of residual element phase error after inter-element 
phase calibration. Fig. 5 shows O, increase in relation to Op. 
Taking into account an overall radiometric error of 1.5K, if 


0.3K is assigned to phase errors, it must be constrained to 0.9° 
residual element phase error. 
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Fig. 5 radiometric error due to element phase error, 

The residual phase error is largely depending on the length of 
averaging time slice for inter-element phase calibration. 
However, the time length cannot be arbitrarily extended due to 
the changing of objective image and instrument stability. After 
one hundred times of repeating experiments, the relationship 


between O,, and averaging time is shown as Fig. 6. 0.9° 


residual element phase error corresponds to 5min averaging 
time. 


0.2 5 10 15 20 25 
averaging time /minute 
Fig. 6 residual element phase error after inter-element phase 
calibration 
From Fig. 5, the sensitivity to phase errors has been found to 
be 0.33 K/deg . However, it is know that the radiometric error 


is mainly contributed by phase errors in the shorter baselines 
due to the fact that it contains a larger amount of signal 
energy''”!, After this inter-element phase calibration, the 
residual baseline phase error is shown in Fig. 7. It should be 
noticed that, the residual phase error is much better for short 
baselines because of the deliberate consideration in the final 
redundancy equations, which is important to reduce the 
radiometric error. 


201703.00288v1 


chinaXiv 


45 


- averaging time: 12.3min “> 
4) yr 
D ome: </mn | 8 
© 3.5} 4 th 
R . fi : 
S 3) g 
E 
o 
@ 2.5} 
a 
£ 
a 2) 
- AF 

È 15r a a 
o T+ 
S 
S 1 


PRE ON 


<x eT 
2 NS fe 


200 300 400 500 600 
baseline length /wavelength 
Fig. 7 residual baseline phase error 
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B. Field Image Calibration 


A proof-of-concept demonstrator with 28 antenna/receiver 
units has already been achieved for GIMS advanced research!" 
as shown in Fig.8. It has the threshold system performance 
requirements proposed by application study, in which, the 
image refresh period is 5mins, and the radiometric accuracy is 
1.5K. 

Imaging tests had been conducted for some angularly 
small targets surrounded by cold sky at a far distance, which is 
beneficial to suppress aliasing effects. A chosen target is shown 
in Fig 9. 

The reconstructed image without phase calibration is shown 
in Fig.10(a). Improved imaging results have been reported in [6] 
by using an external point noise source for phase calibration. 
Here we reprocess the raw data by inter-element phase 
calibration approach. The reconstructed image is demonstrated 
in Fig.10(b), compared with Fig. 10(a), the imaging is improved 


apparently. 
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Fig.8 GIMS proof-of-concept demonstrator Fig.9 optical 


imaging of target 
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Fig.10 (a) raw retrieved image (b) self-calibrated image 


8 
ge 8 @ 


agp. 


V. CONCLUSION 


The traditional phase calibration method of noise injection 
suffers from system complexity and excessive requirements on 
the calibration of noise injection network itself. In this paper, a 


new inter-element phase self-calibration method for rotation 
interferometric radiometers has been introduced, which uses 
the observation scene as reference without needing for extra 
hardware structure and auxiliary point source. 

In the inter-element phase calibration approach, running 
average scheme is proposed to reduce the effect of noise and 
phase ambiguity is solved by preset-phase of very short 
baseline visibility. Simulation and experiments have been 
demonstrated that this method can effectively reduce system 
phase error. 

As future work, some parameter in this new method 
should be optimized for GIMS and calibration precision should 
be further quantificational calculated for the evaluation of 
sensitivity. 
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